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Abstract 

We study, by numerical simulations on a lattice, the dependence of the gauge-invariant 
two-point field-strength correlators in QCD on the path used to perform the color parallel 
transport between the two points. 
(PACS code: 12.38. Gc) 



* Partially supported by MIUR (Italian Ministry of the University and of Scientific and Technological 
Research) and by the INTAS contract 00-0110. 



1. Introduction 



The gauge-invariant two-point correlators of the field strengths in the QCD vacuum are 
defined as 

V^ ua (x) = g^O\Tv{G^ p (0)S(0,x\C)G V(T (x)S\0,x\C)}\0) , (1.1) 

where G w = T a G a w and T a are the matrices of the algebra of the color group SU(3) in 
the fundamental representation The trace in ( [[ . 1\ ) is taken on the color indices. 

S(0,x\C)=Pexp(ig£ dz^z)^ , (1.2) 

with Ap = T a A a , is the Schwinger phase operator needed to parallel-transport the tensor 
G U cr(x) to the point 0. "P" stands for "path ordering" and C is any path from to x. 

Field-strength correlators play an important role in hadron physics. In the spectrum 
of heavy QQ bound states, they govern the effect of the gluon condensate on the level 
splittings [I, 0> P|- They are the basic quantities in models of stochastic confinement of 
color 0, ||, H and in the description of high-energy hadron scattering |7|, ||, || [K| [II], 0, [13 



In some recent works |14|, jlS] , these correlators have been semi-classically evaluated in 
the single-instanton approximation and in the instanton dilute-gas model, so providing 
useful information about the role of the semiclassical modes in the QCD vacuum. 

The correlators ( |1.1[ ) in principle depend on the choice of the path C. Usually, in 
developing the "Stochastic Vacuum Model" (SVM), the path C appearing in Eqs. (Op 
and (|1.2j) is assumed to be the straight line connecting the points and x. Lattice data 



16, p_4 19 , which are the input for these SVM calculations, also refer to the choice 



of the straight-line parallel transport. With this choice the Schwinger string ( |L2|) reads: 

5(0, x) = P exp (ig J dt x^A^xt)) . (1.3) 

Then, in the Euclidean region, translational, 0(4)- and parity invariance require the 
correlator (|1.1|) to be of the following form @, ||, ||] : 



T> mv(r (x) = (5^ u 5p a - S^Spv) V(x 2 ) +T>x 



dVAx 2 ) 

~\~{x ^Xp5 p a x pX(j6 p V -\- XpX a bp iV x pX u 5 p<j) Q^p* ' (^ -Z ^) 
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where T> and T>\ are invariant functions of x 2 . 

The functions D(x 2 ) and T>i(x 2 ) have been directly determined by numerical simula- 



tions on a lattice in the quenched (i.e., pure-gauge) theory, with gauge group SU(2) fll6fl , 
in the quenched SU(3) theory in the range of physical distances between 0.1 and 1 fm 



I7| , |18f and also in full QCD, i.e., including the effects of dynamical fermions |TI 

No analysis exists of what happens for a different choice of the path C than the straight 
line, except for qualitative statements that the physics should not strongly depend on the 
deformations of the path BL H|. In this paper we compute the gauge-invariant field- 
strength correlators for deformed paths (see Fig. 1), in the quenched SU(3) lattice gauge 
theory, for the benefit of the developments of the SVM. 



2. Computations and results 

In principle, every choice of a path C[o,z] connecting the two points and x in the expression 
( |1.2j ) for the Schwinger string operator will generate a different field-strength correlator, 
that we have denoted in Eq. as V^ ua .(x). 

On the lattice we can define a lattice operator T>jQ^ a , which is proportional to V^ ucr 
in the continuum limit, when the lattice spacing a — > 0. Since the lattice analogue of the 
field strength is the open plaquette II w (n) (the parallel transport along an elementary 
square of the lattice lying on the /zp-plane, starting from the lattice site n and coming 
back to n), Vj$£ a will be defined as |T7|, [TJ, |HJ 



V { 2,L(da) = ^{(TT[U^(n)S(n,n + da\C)n ua (n + da)S\n,n + da\C))} , (2.1) 



where 3? stands for real part and the lattice operator VL va (n) is given by |T7|, [18 

n ua {n) = H» - DU(n) - ^Tr[ITt» - IW(n)] . (2.2) 

As explained in Refs. [[17], [L8|], the inclusion of the operator )fl va (in place of nJ, CT ) on the 
right-hand side of Eq. (|2.2|) ensures that the disconnected part and the singlet part of the 
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correlator are left out, so that we are indeed taking the correlation of two operators with 
the exchange of the quantum numbers of a color octet. 

The lattice site n+da is the site at a distance d lattice spacings from n, in the direction 
of a coordinate axis. The Schwinger string S(n, n + da\C) is a parallel transport connecting 
the point n with the point n + da along the path C. 

We have explored the dependence on the path by measuring the correlators for the 
various paths shown in Fig. 1. Fig. la shows the straight-line path [corresponding to the 
straight-line Schwinger operator in Eq. (|1.3|)], which, as explained in the introduction, 



has already been widely analyzed. Fig. lb shows our second choice: the transverse size 
5, in units of the lattice spacing a, parametrizes the deviation from the straight-line case. 
We shall label by b$ the paths in Fig. lb, for a given value of the transverse size 5. If 
an average is made over the orientation of the staple in the plane orthogonal to d, the 
correlator will have the same symmetries (translation, 0(4) and parity) as the one in Fig. 
la and the parametrization ( |1.4|) will still apply, with, of course, different functions X>( bf ) 
and Vf s) . The same considerations also apply to the collection cs of paths shown in Fig. 
lc, for a given value of the transverse size 5, if we average over the orientation of the 
plane containing the path. We have chosen the inversion point p in Fig. lc just in the 
middle of the line (0, x), i.e., p = x/2, in order to preserve parity invariance. For other 
choices, i.e., p = ax with a G [0, 1] and a ^ 1/2, one should average between the path 
with p = ax and the path with p = (1 — a)x in order to preserve parity. 

In both cases, we define |T7[ the two independent functions v\F\x 2 ) and T>^\x 2 ) as 



follows. We go to a reference frame in which x M is parallel to one of the coordinate axes, 
say n = 0. Then 

Vf = \ ± vfjl J (x)=V^+V ( f\ (2.3) 

6 i<j=l 

where = V^ c \x 2 ) and vf ] = vf ] (x 2 ) are the two invariant functions entering in the 
parametrization ( |1.4|) for the given collection of paths C. 

In the naive continuum limit (a — > 0) 

&3Z*fa) a Z a^SU**) + °( a6 ) ■ ( 2 - 4 ) 
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Making use of the definition (|2.3| ) we can also write, in the same limit, 

vf )L (d 2 a 2 ) ~ a 4 vf\d 2 a 2 ) + 0(a 6 ) , 

II a^O II 

V { f )L (d 2 a 2 ) ~ a 4 £>f (dV) + G(a 6 ) . (2.5) 

a— >0 

In order to remove the lattice artefacts, i.e., the terms 0(a 6 ) in Eqs. (|2.4j) and (|2.5|) , we 
shall make use of the cooling technique, described in previous papers (see Refs. |20], |2T|| 
and |I7|, [IE], Cooling is a local procedure, which affects correlations at distances 



that grow as the square root of the number of cooling steps, as in a diffusion process. 
If the distance at which we observe the correlation is sufficiently large, we then expect 
that lattice artefacts are frozen by cooling long before the correlation is affected: this will 
produce a plateau in the dependence on the cooling step. Our data are the values of the 
correlators at the plateau; the error is the typical statistical error at the plateau, plus a 
systematic error which is estimated as the difference between neighbouring points at the 
plateauM 

We have measured the correlations on a 16 4 lattice at distances ranging from 3 to 8 
lattice spacings and at f3 — 6. The value of the lattice spacing a in physical units can be 
extracted from the value of the string tension [p2| [23| and it is a ~ 0.1 fm at our value 
of (3. At this value of j3 scaling has been already tested for the case of the straight-line 
paths in Refs. [|Tj], [Tj|. The data of Refs. [|17| , |18f came from several different values 



of (3 (including also (3 — 6) and two lattice sizes, 16 4 and 32 4 . They showed no visible 
dependence neither on the ultraviolet, nor on the infrared scale. We therefore do not 
expect significant scaling violations. 

The results are shown in Figs. 2, 3, 4 and 5, versus the distance d in units of the 
lattice spacing a. The lines are drawn as an eye-guide. Figs. 2 and 3 refer to the paths 
in Fig. lb, for various values of the transverse size 5 of the staple. Figs. 4 and 5 give the 
analogous results for the paths shown in Fig. lc. 

The very result of this paper is the unexpectedly strong dependence of the correlators 
on the shape of the path. The correlator with the straight-line path ( |1.3| ) (see Fig. la), 
has the largest signal (for every distance d), compared with the two other choices for the 
path. Every deformation from the straight-line path seems to produce a sharp decrease 
of the value of the correlator. What seems to be path-independent is the slope of the 

tin Figs. 2, 3, 4 and 5 we have plotted only the points corresponding to a clear plateau in the cooling 
process. For some points at large distances, our cooling proved to be not long enough to reach the plateau. 
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curves, which look parallel to one another in the linear-log plots of Figs. 2, 3, 4 and 5. 
This means that the correlation length A^ defined in Refs. [|17], [18], |19[ is approximately 
path-independent, at least for the classes of paths that we have considered. 

In Figs. 2, 3, 4 and 5 the best fit to the data for straight-line Schwinger string of 
Ref. fl8] , Eqs. (2.10) and (2.11), is displayed for comparison (see also Ref. [2Hfl , where a 
critical comparison among different best fits is performed). 

In order to quantitatively test the independence of the correlation length Xa on the 
path C, we have tried a best fit to the data for T>± — T> + T>i, at intermediate distances 
5 < d < 8, i.e., about 0.5 fm < x < 0.8 fm), with the same function used in Ref. | 18| , Eq. 
(2.10): 



Pi 



[x 



A ± e 



-\x\/x A , a±_ 



\x\/\ a 



(2.6) 



4.9 



Keeping A^ fixed to the value of Ref. Eq. (2.11), i.e., Xa = 1/182Al, where 
MeV is the lattice scale [2^, |23" |, a good fit results with a reasonable x 2 /^d.o.f.- The path- 
dependence of the correlator reflects in a rather strong path-dependence of the coefficient 
A± of the exponential term. In the range of distances chosen, the second term of Eq. 
(|2.6| ) is compatible with zero, within the errors. The coefficient A± changes up to a factor 
of four, going from 5 = to 5 = 6 for the collection of paths b$ (Fig. lb); and it changes 
up to a factor of two, going from 5 = to 5 = 6 for the collection of paths cs (Fig. lc). 

In Ref. [0 Eq. ( pHf ) was considered as a split-point regulator of the gluon condensate, 
which was extracted from the coefficient Aj_ of the exponential term (see also Ref. P^]). 
It is not clear to us how this determination depends on the shape of the path. 
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FIGURE CAPTIONS 



Fig. 1. The different paths for the Schwinger string S that we have considered for our 
analysis. 

Fig. 2. The function a A T>^^ versus the distance d in units of the lattice spacing a. Circles 
correspond to the straight-line path of Fig. la; triangles to Fig. lb with 5 = 1; 
squares to Fig. lb with 5 = 4; diamonds to Fig. lb with 5 = 6. Lines are drawn 
as eye-guides. The thick continuum line has been obtained using the parameters of 



the best fit obtained in Ref. P], Eqs. (2.10) and (2.11) 



Fig. 3. The function a 4 V^ versus the distance d in units of the lattice spacing a. The 
symbols are the same as in Fig. 2. 

Fig. 4. The function a 4 V^ versus the distance d in units of the lattice spacing a. The 
symbols are the same as in Fig. 2, except that they now refer to the path of Fig. 
Ic. 

Fig. 5. The function a 4 V^ versus the distance d in units of the lattice spacing a. The 
symbols are the same as in Fig. 2, except that they now refer to the path of Fig. 
Ic. 
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Figure 1 
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